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Abstract 

Using  a  middle  atmospheric  general  circulation  model,  we  find  that  zonally  asymmetric 
ozone  (ZAO)  profoundly  changes  the  morphology  of  the  Northern  Hemisphere  planetary 
waveguide  (PWG).  ZAO  causes  the  PWG  to  contract  meridionally  and  expand  vertically, 
with  a  significant  increase  in  wave  propagation.  Consequently,  there  is  a  significant 
increase  in  the  upward  flux  of  wave  activity  from  the  troposphere  and  lower  stratosphere 
into  the  interior  of  the  stratosphere  and  lower  mesosphere.  The  ZAO-induced  changes  in 
the  PWG  increase  the  Eliassen-Palm  flux  divergence,  causing  a  warmer  and  weaker 
stratospheric  polar  vortex.  The  ability  of  ZAO  to  alter  the  flux  of  planetary  wave  activity 
into  the  polar  vortex  has  significant  implications  for  accurately  modeling  wave-driven 
phenomena  in  the  middle  atmosphere,  including  sudden  stratospheric  warmings,  11 -year 
solar  cycle-modulated  wave  activity,  and  the  Brewer-Dobson  circulation. 

1.  Introduction 

Understanding  the  physics  that  controls  the  planetary  waveguide  (PWG)  is  not 
fully  understood;  yet  such  understanding  is  vital  to  producing  reliable  assessments  of 
human-caused  impacts  on  key  features  of  the  climate  system,  including  the  upward  flux 
of  planetary  wave  activity,  the  frequency  of  stratospheric  sudden  warmings,  and  the 
Brewer-Dobson  circulation.  Here  we  demonstrate  that  the  physics  associated  with  zonally 
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asymmetric  ozone  (ZAO)  profoundly  changes  the  morphology  of  the  PWG,  a  result  that 
underscores  the  importance  of  accurately  accounting  for  ZAO  in  global  climate  models. 

The  concept  of  a  PWG  was  first  discussed  by  Dickinson  [1968],  who  proposed 
that  planetary  waves  would  tend  to  propagate  upward  through  the  region  of  weak 
westerly  winds  bounded  by  the  North  Pole  and  the  strong  westerly  winds  of  the 
stratospheric  polar  vortex.  Matsuno  [1970]  derived  a  refractive  index  squared  (RI)  for 
planetary  waves  on  a  sphere  and  showed  that  while  the  large  westerly  wind  speeds  within 
the  polar  night  jet  may  impede  wave  propagation,  the  vortex  edge  is  characterized  by 
strong  meridional  gradients  in  planetary  vorticity  that  enhance  wave  propagation.  Thus, 
rather  than  the  polar  vortex  acting  as  the  southern  boundary  of  the  waveguide,  it  is  the 
tropical  zero  wind  line  that  limits  the  equatorward  propagation  of  planetary  waves  [Tung, 
1979].  The  PWG,  therefore,  is  largely  defined  by  the  strength  of  wave  propagation 
(measured  by  the  RI)  and  its  meridional  width  (measured  by  the  location  of  the  zero  wind 
line).  Mechanistic  modeling  studies  have  shown  that  diabatic  effects  due  to  ZAO 
modulate  both  of  these  and  other  wave-driven  features  in  the  tropics  [Echols  and  Nathan, 
1996;  Cordero  and  Nathan,  1998,  2000]  and  extratropics  [Nathan  and  Li,  1991;  Nathan 
et  al.,  1994;  Nathan  and  Cordero,  2007]. 

For  example,  Nathan  and  Cordero  [2007]  derived  an  ozone-modified  refractive 
index  (OMRI)  for  vertically  propagating  planetary  waves  that  describes  how  ZAO  can 
produce  changes  in  planetary  wave  propagation  and  wave  attenuation  (damping). 
Specifically,  the  OMRI  makes  clear  how  ZAO,  which  is  a  function  of  the  vertical  and 
horizontal  gradients  of  background  ozone,  modulates  planetary  wave  drag.  Depending  on 
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the  background  ozone  gradients,  planetary  wave  drag  may  either  increase  or  decrease  due 
to  ZAO-induced  changes  in  wave  propagation  and  wave  damping. 

Evidence  supporting  the  importance  of  ZAO  in  the  stratosphere  is  not,  however, 
limited  to  mechanistic  studies.  ZAO  has  also  been  shown  to  be  important  in  altering  the 
climate  of  general  circulation  models  [e.g.,  Sassi  et  al,  2005;  Gabriel  et  al.,  2007;  Crook 
et  al,  2008;  Waugh  et  al,  2009;  Gillett  et  al,  2009;  McCormack  et  al,  2011].  For 
example,  McCormack  et  al  [2011]  used  a  general  circulation  model  with  fully  interactive 
ozone  and  dynamics  to  demonstrate  that  ZAO  produces  a  warmer  and  weaker  Northern 
Hemisphere  stratospheric  polar  vortex,  along  with  an  increased  incidence  of  SSWs. 

Neither  the  mechanistic  studies  nor  the  general  circulation  model  studies  cited 
above  have  addressed  the  effects  of  ZAO  on  the  PWG.  Owing  to  the  importance  of 
accurately  representing  the  PWG  in  global  climate  models,  the  goal  of  this  paper  is  to 
examine  how  ZAO  alters  the  morphology  of  the  PWG  of  the  Northern  Hemisphere 
during  winter.  In  the  following  section  we  introduce  the  model  and  diagnostics  that  we 
use  to  measure  changes  in  the  morphology  of  the  PWG.  We  then  present  the  results  and 
follow  with  a  discussion  of  the  implications  for  global  climate  modeling. 

2,  Model  Description  and  Diagnostics 

We  use  the  same  model  and  experimental  setup  as  in  McCormack  et  al  [2011]. 
The  effects  of  ZAO  on  planetary  wave  propagation  during  Northern  Hemisphere  winter 
are  examined  using  the  NOGAPS-ALPHA  global  spectral  general  circulation  model  [see, 
e.g.,  Eckermann  et  al,  2009;  McCormack  et  al,  2009  and  references  therein],  which 
extends  from  the  surface  to  ~90  km  in  height  and  utilizes  68  hybrid  (o-p)  vertical  levels 
with  triangular  truncation  at  wavenumber  79.  NOGAPS-ALPHA  includes  prognostic 
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equations  for  O3  and  H2O  calculated  using  the  photochemical  parameterizations  of 
McCormack  et  al.  [2006,  2008].  Shortwave  heating  and  longwave  cooling  rates  are 
computed  using  the  predicted  profiles  of  O3  and  H2O,  and  a  fixed  profile  of  CCK  The 
lower  boundary  condition  is  computed  from  observed  12-hourly  sea  surface  temperature 
and  surface  ice  distributions. 

In  total,  we  ran  fifteen  pairs  of  simulations  spanning  early  December  to  late 
March.  Each  pair  of  model  simulations  was  initialized  using  identical  profiles  of  wind, 
temperature,  and  chemical  constituents  acquired  from  the  NOGAPS-ALPHA  data 
assimilation  system  [ Eckermann  et  al.,  2009].  The  data  we  used  for  each  pair  of 
simulations  was  staggered  in  time  in  order  to  generate  an  ensemble  of  fifteen  independent 
pairs  of  model  simulations  (for  details  see  McCormack  et  al.,  2011).  For  each  pair  of 
simulations,  one  of  the  simulations  utilized  fully  prognostic  ozone  in  the  radiative  heating 
and  cooling  calculations  (designated  3D03),  and  the  other  simulation  utilized  zonal-mean 
ozone  values  to  evaluate  the  radiative  heating  and  cooling  rates  (designated  ZM03). 
Zonal-mean  ozone  values  were  calculated  from  the  prognostic  ozone  field.  We  isolate  the 
effect  of  ZAO  on  the  circulation  by  taking  the  difference  between  3D03  and  ZM03 
simulations  for  any  given  variable.  Statistical  significance  between  model  variables  was 
then  assessed  using  a  Student’s  T-test  at  the  95%  confidence  interval.  As  stated  in 
McCormack  et  al.  [2011],  the  3D03  model  runs  produced  four  SSWs,  while  the  ZM03 
runs  produced  only  one  SSW.  Statistical  significance  calculations  were  also  carried  out 
with  all  of  the  runs  containing  SSWs  removed;  however,  doing  so  did  not  change  the 
qualitative  nature  of  our  results. 
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The  morphology  of  the  PWG  is  evaluated  in  terms  of  two  properties:  i)  its  shape 
(meridional  width  and  vertical  extent);  and  ii)  the  strength  of  propagation  within  the 
guide  itself.  To  diagnose  the  morphology  of  the  PWG  and  to  distinguish  regions  of  wave 
propagation  versus  evanescence,  we  employ  the  spherical  fonn  of  the  quasigeostrophic 
refractive  index  squared  (RI)  [Andrews  et  al.,  1987]: 


(1) 


where 


2Q  cos  (f)  if  {u  cos  (/>\ 
a  <7 2  I,  cos  (f) 


f  uz 

1 -  Ps~i 

Ps  l  N2 


is  the  meridional  potential  vorticity  gradient;  <f>  is  latitude,  z  is  height,  k  is  the  zonal 
wavenumber,  N  is  the  buoyancy  frequency,  /  is  the  Coriolis  parameter,  H  is  the  mean 
scale  height  [=7  km],  ps  is  the  sea-level  reference  density,  a  is  the  radius  of  the  Earth, 

u  is  the  zonal  mean  zonal  wind,  Q  is  the  Earth’s  rotation  frequency,  and  subscripts 
denote  derivatives  with  respect  to  the  given  variable. 

Planetary  waves  propagate  within  regions  where  n2k  >  0  and  are  evanescent  in 

regions  where  nl<  0 .  Because  \n\  |  — »  oo  as  u  — »  0 ,  interpretation  of  the  RI  near  the  zero 
wind  line  becomes  problematic.  What  is  important  here  is  that  near  the  zero  wind  line  the 
sign  of  u  largely  controls  the  sign  of  n\ .  We  therefore  use  the  location  of  the  zero  wind 
line  to  approximate  the  southern  boundary  (and  meridional  width)  of  the  PWG.  Our  use 
of  the  zero  wind  line  as  a  measure  of  the  southern  boundary  of  the  PWG  is  consistent 
with  previous  studies  [e.g.,  Holton  and  Tan,  1982].  The  vertical  extent  of  the  PWG  and 
the  strength  of  wave  propagation  within  the  waveguide  are  evaluated  using  the  RI. 
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Moreover,  if  waves  are  assumed  to  be  WKB  in  nature,  then  the  local  direction  and 
magnitude  of  the  vectors  of  the  Eliassen-Pahn  flux  (EP-flux)  are  also  related  to  the  RI 
(i.e.,  regions  of  larger  RI  are  associated  with  larger  EP-flux  vectors)  [ Edmon  et  al.,  1980; 
Butchart  et  al.,  1982].  We  plot  EP-flux  vectors  at  100  hPa  to  assess  how  changes  in  wave 
propagation  modulate  the  upward  flux  of  wave  activity  into  the  middle  and  upper 
stratosphere.  When  plotting  the  EP-flux  vectors,  we  adopt  the  scaling  of  Butchart  et  al. 
[1982]  so  that  the  EP-flux  and  its  divergence  are  graphically  consistent  when  plotted  in 
Cartesian  coordinates.  We  measure  the  net  effect  of  changes  in  the  zonal-mean 
circulation  due  to  changes  in  the  EP-flux  divergence  by  the  deceleration  tenn  defined  in 
Palmer  [1981]  as 

exp(zAf/)vifj  (3) 

/?sacos^ 

where  V  •  F  is  the  EP-flux  divergence  and  other  variables  are  defined  above. 

3.  Results 

Figures  la  and  lb  show  the  difference  in  the  ensemble  average  zonal-mean  zonal 
wind  (3D03  minus  ZM03)  for  January  and  February;  wind  differences  that  are 
statistically  significant  above  the  95%  confidence  interval  are  shaded.  Because  the  largest 
change  in  wind  speed  occurs  in  February  in  the  Northern  Hemisphere  (NH),  we  focus  our 
attention  on  this  month  and  region  for  the  remainder  of  the  paper. 

As  Figure  1  shows,  ZAO  produces  statistically  significant  decreases  in  the  zonal 
mean  westerly  flow  of  5-11  ms'1  in  the  equatorial  upper  stratosphere  and  middle 
mesosphere,  as  well  as  decreases  of  3-18  ms'1  throughout  the  NH  extratropical 
stratosphere  and  mesosphere.  The  difference  in  the  NH  zonal-mean  wind  produces  two 
notable  results.  First,  the  decrease  in  the  wind  speed  in  the  tropical  upper  stratosphere  and 
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lower  mesosphere  causes  a  northward  migration  of  the  tropical  easterlies  and  a  northward 
shift  of  the  zero  wind  line;  the  northward  shift  in  the  zero  wind  line  varies  from  1°  to  8° 
latitude  between  30  km  and  65  km  in  height  (Figure  2).  Second,  the  decrease  in  wind 
speed  within  the  extratropics  and  polar  region  represents  a  weaker  polar  vortex  and  a 
wanner  stratosphere  (see  also  McCormack  et  al.  [2011],  their  Figure  2).  These  two 
changes  in  wind  structure  combine  to  modify  the  morphology  of  the  PWG  in  distinct 
ways;  the  change  in  the  location  of  the  zero  wind  line  narrows  the  width  of  the 
waveguide,  while  the  change  in  the  extratropical  and  polar  winds  expands  the  PWG  in 
the  vertical  and  strengthens  wave  propagation  within  the  waveguide. 

The  meridional  contraction  of  the  PWG  and  the  warmer,  more  disturbed  polar 
vortex  found  in  our  results  is  in  qualitative  agreement  with  recent  studies  suggesting  a 
connection  between  northward  shifts  in  the  upper  stratospheric  zero  wind  line  and  an 
increase  in  the  frequency  of  SSWs.  For  example,  Gray  [2003]  found  that  in  a  similar 
fashion  to  the  so-called  ‘Holton-Tan’  mechanism  in  the  lower  stratosphere  [ Holton  and 
Tan,  1982],  a  northward  shift  in  the  upper  stratospheric  zero  wind  line  produces  a 
narrower  waveguide  that  channels  more  wave  activity  poleward.  The  increase  in  wave 
activity  induces  more  wave  drag  on  the  zonal-mean  circulation  and  creates  a  warmer 
extratropical  and  polar  stratosphere  and  an  increased  incidence  of  SSWs.  Indeed,  as  our 
EP-flux  diagnostics  will  show,  the  combination  of  a  narrower  PWG  and  an  increase  in 
wave  propagation  (discussed  next)  combine  to  produce  an  increase  in  wave  activity  and 
the  weaker  polar  vortex  shown  in  Figure  lb. 

Figures  3a  and  3b  compare  the  contours  of  the  February  ensemble  average  RI  (1) 
for  ZM03  and  3D03  for  stationary  planetary  wave  (s=l);  similar  results  were  obtained 
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for  planetary  wave  (s=2).  The  RI  scale  extends  from  red  to  blue,  which  corresponds  with 
large  to  small  RI  values,  respectively;  regions  where  planetary  waves  are  evanescent 
(RI<0)  are  denoted  by  white  space.  Figure  3  shows  that  ZAO  increases  planetary  wave 
propagation  in  two  important  ways.  First,  the  vertical  extent  of  the  PWG  increases,  which 
allows  the  waves  to  propagate  10  to  15  km  higher  within  the  mesosphere  between  50°- 
70°.  Second,  ZAO  causes  a  notable  increase  in  the  magnitude  of  the  RI  throughout  the 
entire  stratosphere.  The  latter  point  can  be  seen  more  easily  by  examining  the  difference 
in  RI  (3D03  minus  ZM03)  shown  in  Figure  3c.  The  change  in  the  RI  (Figure  3c)  shows 
that  ZAO  produces  a  significant  increase  in  wave  propagation  from  ~15  km  to  80  km  in 
height  between  30°-70°  north  latitude.  To  better  understand  the  implications  of  the 
expansion  of  the  PWG  in  the  vertical  and  the  increase  in  wave  propagation  within  the 
waveguide,  we  next  consider  the  EP-flux  and  its  divergence. 

Figures  4a-c  show  contours  of  the  EP-flux  divergence  in  February  for  ZM03, 
3D03,  and  the  difference  3D03  minus  ZM03.  In  addition,  each  of  the  EP-flux 
divergence  contour  plots  is  overlaid  with  the  corresponding  EP-flux  vectors  depicting  the 
direction  and  relative  magnitude  of  wave  propagation  at  100  hPa.  The  EP-flux  vectors 
displayed  in  Figure  4  clearly  show  that  ZAO  causes  an  increase  in  the  upward  flux  of 
wave  activity  between  ~30°-75°;  EP-flux  vectors  plotted  at  levels  throughout  the 
stratosphere  show  a  similar  increase  due  to  ZAO.  As  expected,  the  increase  in  upward 
wave  activity  corresponds  remarkably  well  with  the  increase  in  RI  shown  in  Figure  3c. 
Moreover,  the  region  of  increased  upward  wave  activity  depicted  by  the  EP-flux  vectors 
is  also  associated  with  a  larger  and  stronger  region  of  EP-flux  divergence  below  30  km 
(Figure  4c).  These  results  indicate  that  changes  in  wave  propagation  due  to  ZAO  within 
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the  lower  stratosphere  play  an  important  role  in  modulating  the  amount  of  wave  activity 
that  is  allowed  to  propagate  from  the  troposphere  into  the  interior  of  the  stratosphere  and 
lower  mesosphere.  The  importance  of  the  lower  stratospheric  changes  in  EP-flux  and  its 
divergence  is  further  corroborated  by  considering  the  change  in  the  zonal-mean 
deceleration  (3).  Figures  5a-c  show  the  deceleration  of  the  zonal-mean  wind 
corresponding  to  the  EP-flux  divergence  of  Figure  4  for  ZM03,  3D03,  and  the  difference 
3D03  minus  ZM03.  Indeed,  Figure  5c  shows  a  substantial  net  increase  in  the  zonal- 
mean  deceleration,  Aj,  due  to  ZAO.  The  increase  in  Af  correlates  very  well  with  both  the 
expansion  of  the  PWG  in  the  vertical  and  the  increase  in  wave  propagation  within  the 
waveguide.  The  increase  in  Af  also  corresponds  very  well  with  the  decrease  in  the  NH 
stratospheric  and  mesospheric  zonal-mean  wind  profile  shown  in  Figure  lb. 

4.  Conclusions  and  discussion 

We  investigated  the  effect  of  zonally  asymmetric  ozone  (ZAO)  on  the  zonal -mean 
circulation  and  planetary  waveguide  (PWG)  of  the  Northern  Hemisphere  winter 
stratosphere  and  lower  mesosphere  using  an  ensemble  of  global  circulation  model  (GCM) 
simulations.  We  demonstrated  that  ZAO  profoundly  alters  the  morphology  of  the  PWG  in 
two  distinct  ways.  First,  the  PWG  contracts  meridionally  and  expands  vertically  within 
the  upper  stratosphere  and  lower  mesosphere.  Second,  the  magnitude  of  wave 
propagation  increases  throughout  the  interior  of  the  stratospheric  waveguide  between 
~30°-75°.  In  combination,  the  changes  in  the  shape  and  strength  of  the  waveguide  are 
associated  with  an  increase  in  the  deceleration  of  the  zonal-mean  circulation,  a  weaker 
and  wanner  stratospheric  polar  vortex,  and  an  increased  incidence  of  SSWs. 
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The  modeled  changes  in  the  PWG  have  important  implications  for  modeling 
SSWs,  variations  in  stratospheric  dynamics  associated  with  the  11-year  solar  cycle,  and 
the  Brewer-Dobson  circulation.  For  example,  Charlton  et  al.  [2007]  investigated  the 
ability  of  a  series  of  GCMs  -  which  did  not  consider  the  effects  of  ZAO  on  the  PWG  -  to 
generate  an  accurate  number  of  SSWs;  the  results  showed  that  most  GCMs  failed  to 
produce  enough  SSWs  when  compared  to  observations.  In  the  current  GCM  study  - 
where  we  explicitly  consider  the  effects  of  ZAO  on  the  PWG  -  the  model  runs  with  ZAO 
had  a  noticeably  larger  EP-flux  at  100  hPa  and  produced  a  larger  number  of  SSWs 
compared  to  similar  runs  without  ZAO.  This  increase  in  the  EP-flux  correlates  very  well 
with  the  increase  in  the  RI  in  the  lower  stratosphere,  indicating  that  changes  in  the  PWG 
associated  with  ZAO  may  play  an  important  role  in  modulating  the  flux  of  wave  activity 
entering  the  stratosphere  from  the  troposphere,  which  is  known  to  drive  SSWs. 

SSWs  have  also  been  correlated  with  northward  shift  of  tropical  easterlies  and  the 
zero  wind  line  within  the  upper  stratosphere  and  lower  mesosphere  [Gray,  2003;  Vineeth 
et  al.,  2010].  Gray  et  al.  [2003]  suggested  that  because  vertically  propagating  planetary 
waves  have  vertical  wavelengths  comparable  to  the  distance  between  the  tropopause  and 
the  upper  stratosphere  during  mid-  to  late-winter,  it  is  possible  that  shifts  of  the  zero  wind 
line  in  the  upper  stratosphere  help  to  guide  waves  poleward  in  much  the  same  manner  as 
the  traditional  Holton-Tan  mechanism  in  the  lower  stratosphere.  Although  our  results 
have  not  tested  the  actual  mechanics  of  this  theory,  it  is  nonetheless  interesting  to  note  the 
correlation  between  the  northward  shift  in  the  zero  wind  line  and  the  increase  in 
frequency  of  SSWs  in  our  model  runs  that  include  ZAO. 
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In  addition  to  SSWs,  two  aspects  of  the  PWG  discussed  in  this  study  -  the  width 
and  strength  -  may  provide  a  means  for  ZAO  to  communicate  and  amplify  the  signal 
from  the  1 1-year  solar  cycle.  The  relationship  between  the  width  of  the  PWG  and  solar- 
modulated  ZAO  has  been  shown  by  Cordero  and  Nathan  [2005]  to  be  an  important 
pathway  for  communicating  the  solar  signal  to  the  QBO  (see  their  Figure  1).  Because  our 
results  show  that  ZAO  plays  a  role  in  modulating  the  upper  stratospheric  tropical 
easterlies,  it  is  therefore  likely  that  in  a  similar  manner  to  the  Cordero  and  Nathan  solar- 
QBO  mechanism,  ZAO  may  also  play  a  role  in  modulating  the  semi-annual  oscillation 
(SAO)  in  the  tropical  upper  stratosphere.  Because  the  SAO  plays  a  dominant  role  in  the 
location  and  timing  of  shifts  in  the  upper  stratospheric  zero  wind  line,  modulation  of  the 
SAO  by  ZAO  may  provide  an  important  pathway  for  solar  modulation  of  the  width  of  the 
PWG.  Indeed,  this  is  a  topic  that  deserves  further  attention  given  the  parallel  timing  of 
the  peak  in  wind  modulation  from  ZAO  in  our  results  and  the  robust  solar  signal 
observed  during  mid-  to  late-winter  in  the  tropical  upper  stratosphere  [e.g.,  Labitzke  and 
van  Loon,  1988;  Naito  and  Hirota,  1997].  In  addition  to  solar  modulation  of  the  width  of 
the  waveguide,  ZAO  also  provides  a  way  to  change  the  strength  of  wave  propagation 
within  the  waveguide  itself.  This  scenario  is  supported  by  recent  work  demonstrating  the 
link  between  the  1 1  -year  solar  cycle  and  changes  in  wave  propagation  and  planetary 
wave  drag  in  the  extratropical  stratosphere  [ Nathan  et  al,  20 11], 

Finally,  we  note  that  recent  modeling  studies  have  projected  that  increasing 
greenhouse  gas  concentrations  will  affect  stratospheric  ozone  abundance  [Eyring  et  al., 
2007]  and  produce  increases  in  wave  activity  that  will  accelerate  the  Brewer-Dobson 
circulation  [Garcia  and  Randel,  2008].  Because  ZAO  is  intimately  connected  to  ozone 
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abundance  and  wave  activity,  ZAO  may  provide  an  important  additional  pathway  for 
communicating  greenhouse  gas-induced  changes  in  stratospheric  ozone  to  the  Brewer- 
Dobson  circulation.  We  are  currently  investigating  these  possibilities. 
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(a)  Change  in  Zonal-Mean  Zonal  Wind  for  January 
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(b)  Change  in  Zonal-Mean  Zonal  Wind  for  February 
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Figure  1.  Monihly  average  zonal-mean  zonal  wind  difference  (AU  =  3D03  minus 
ZM03)  for:  (a)  January  and  (b)  February.  Contours  intervals  noted  in  the  Figure  are  in 
units  of  ms'1.  Statistically  significant  wind  changes  at  the  95%  confidence  interval  are 
shaded  in  gray. 
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Figure  2.  Location  of  Northern  Hemisphere  zero  wind  line  for  3D03  (solid  line)  and 
ZM03  (dashed  line). 
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(a)  Rl  February  ZM  Ozone  (s=l)  (c)  a  Rl  February 
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Figure  3.  MonLhly  averaged  refractive  index  squared  (Rl)  in  February  for:  (a)  ZM03  for 
planetary  wave  s=l;  (b)  3D03  for  planetary  wave  s=l;  and  (c)  difference  in  Rl  (ARI  = 
3D03  minus  ZM03).  Regions  of  wave  evanescence  (RI<0)  are  depicted  by  regions 
without  shading  (white).  In  all  plots,  the  Rl  has  been  multiplied  by  a  ;  the  contour 
interval  is  4. 
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(a)  February  ZM  EP  Flux  Oiv  (c)  February  a  EP  Flux  Div 


Latitude  Latitude 


Figure  4.  Monthly  averaged  EP-flux  divergence  in  February  for:  (a)  ZM03;  (b)  3D03; 
and  (c)  difference  in  EP-flux  divergence  (3D03  minus  ZM03).  Arrows  represent  EP- 
vectors  plotted  at  the  100  hPa  height  level  (~16  km).  The  contour  interval  is  5x10'  kg  s'  . 
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Figure  5.  Contours  of  monthly  averaged  zonal-mean  deceleration  (A/)  in  February  for:  (a) 
ZM03;  (b)  3D03;  and  (c)  the  difference  (3D03  minus  ZM03).  The  contour  interval  is 
4xl05  ms'2. 
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